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ABSTRACT 

Recent observational and theoretical progress has favored merging and helium-accreting sub- 
Chandrasekhar mass white dwarfs in the double-degenerate and the double-detonation channels, 
respectively, as the most promising progenitors of normal Type la supernovae (SNe la). Thus the 
fate of rapidly-accreting Chandrasekhar mass white dwarfs in the single-degenerate channel remains 
more mysterious then ever. In this paper, we clarify the nature of ignition in Chandrasekhar-mass 
single-degenerate SNe la by analytically deriving the existence of a characteristic length scale which 
establishes a transition from central ignitions to buoyancy-driven ignitions. Using this criterion, com¬ 
bined with data from three-dimensional simulations of convection and ignition, we demonstrate that 
the overwhelming majority of ignition events within Chandrasekhar-mass white dwarfs in the single¬ 
degenerate channel are buoyancy-driven, and consequently lack a vigorous deflagration phase. We 
thus infer that single-degenerate SNe la are generally expected to lead to overluminous I991T-like 
SNe la events. We establish that the rates predicted from both the population of supersoft X-ray 
sources and binary population synthesis models of the single-degenerate channel are broadly consis¬ 
tent with the observed rates of overluminous SNe la, and suggest that the population of supersoft 
X-ray sources are the dominant stellar progenitors of SNe 1991T-like events. We further demonstrate 
that the single-degenerate channel contribution to the normal and failed 2002cx-like rates is not likely 
to exceed 1% of the total SNe la rate. We conclude with a range of observational tests of overluminous 
SNe la which will either support or strongly constrain the single-degenerate scenario. 

Subject headings: supernovae: general — supernovae: individual (1991T, 2002cx) — hydrodynamics 
— white dwarfs 


1. INTRODUCTION 

Type la supernovae (SNe la) play an important role 
as standardizable candles for cosmology (Schmidt et al. 
1998; Riess et al. 1998; Perlmutter et al. 1999). Yet, the 
nature of their stellar progenitors remains elusive (Maoz 
et al. 2014). Several leading models have been advanced. 
In the classic single-degenerate channel, a hydrogen-rich 
main sequence or red giant star accretes matter onto a 
white dwarf companion, until the companion nears or ex¬ 
ceeds the Chandrasekhar mass and explodes (Whelan & 
Iben 1973). In contrast, the double-degenerate channel 
posits that double C/0 white dwarf systems, brought 
together by gravitational waves, merge and lead to SNe 
la. Additionally, WDs accreting helium from either a 
non-degenerate or WD companion may also give rise to 
a SN la through the double-detonation channel (Nomoto 
1982). Recent observational evidence has begun to favor 
double-degenerates, or possibly double-detonations, as 
the origin for the majority of SNe la. On the other hand, 
additional evidence - ranging from bimodality in spec¬ 
tral features (Benetti et al. 2005), circumstellar material 
(Dilday et al. 2012), and the ejected mass-®®Ni relation 
(Scalzo et al. 2014), as well as stable iron peak-element 
nucleosynthesis (Seitenzahl et al. 2013; Yamaguchi et al. 
2015) - has hinted that multiple channels may contribute 
to the overall SNe la rate. 

Therefore, even if double-degenerate or double- 
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detonation supernovae do account for the majority of 
normal SNe, a number of major questions remain unre¬ 
solved. Specifically, what is the fate of rapidly-accreting 
white dwarfs in the single-degenerate channel once they 
reach the Chandrasekhar mass? If the single-degenerate 
channel does give rise to SNe la, how can we reconcile 
the existence of single-degenerate supernovae with the 
general absence of observational evidence in their sup¬ 
port? Perhaps most crucially, if there are indeed multiple 
channels producing SNe la, what are the specific observa¬ 
tional characteristics which would enable us to separate 
the population of single-degenerate SNe la from the pop¬ 
ulation of SNe la as a whole, and most directly confront 
observations with theory? 

To address these questions, it is important to recall 
that SNa la are standardizable candles, but exhibit a 
range of diverse outcomes, from failed SN 2002cx-like 
events, through normal SNe la, to overluminous SN 
1991T-like SNe la (Filippenko et al. 1992), and super- 
luminous supernovae such as 2003fg (Howell et al. 2006) 
or 2009dc (Yamanaka et al. 2009; Silverman et al. 2011). 
In this paper, we will explore the possible connection of 
single-degenerate SNe la to the population of overlumi¬ 
nous SNe la. Overluminous SNe la account for 4 - 20% 
of all SNe la (Foley et al. 2009; Li et al. 2001), and yield 
a significantly larger amount of ^®Ni than normal SNe la 
(roughly 1 Mq, as opposed to 0.5 — O.SMq). 

The rates and high ®®Ni yields of overluminous SNe la 
are in tension with both double-detonation and double¬ 
degenerate merger models of SNe la, since both fun¬ 
damentally rely upon sub-Chandrasekhar mass WDs. 
The prompt detonation of cold, sub-Chandrasekhar C/0 
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WDs up to 1.15M0 produce only up to .81Mq of ®®Ni, 
which is at the upper end of the range of °®Ni pro¬ 
duced in normal SNe la (Sim et al. 2010). Therefore, 
the production of the amounts of ®®Ni required for over- 
luminous events in either the violent double-degenerate 
or double-detonation scenario requires massive C/0 pri¬ 
maries in excess of I. 2 M 0 , and approaching the Chan¬ 
drasekhar mass, which both observation and modeling 
suggest are very rare in both channels (Badenes & Maoz 
2012; Ruiter et al. 2011, 2013). It may be possible for two 
sub-Chandrasekhar white dwarfs to merge into a near- 
Chandrasekhar mass white dwarf on a viscous evolution¬ 
ary timescale and produce an overluminous SN la (van 
Kerkwijk et al. 2010; Piro et al. 2014), though it is re¬ 
mains unclear precisely whether such mergers will deto¬ 
nate (Shen et al. 2012; Schwab et al. 2012; Ji et al. 2013). 
Furthermore, in the context of double-detonations, mod¬ 
els suggest that even the minimal mass of a thin detonat¬ 
ing shell of helium will in turn detonate the underlying 
C/0 WD core (Fink et al. 2010). Consequently, He- 
accreting massive C/0 WDs may tend to detonate eas¬ 
ily, with little opportunity to accrete the requisite mass 
to explain overluminous SNe la. 

In contrast, overluminous supernovae have long been 
understood to be the natural consequence of the detona¬ 
tion of Chandrasekhar-mass WDs (Arnett 1969). How¬ 
ever, until recently, conventional lore has held that the 
near-uniformity of normal SNe la strongly favored the 
single-degenerate channel as the progenitors of normal 
SNe la. For decades, theorists have worked to rec¬ 
oncile the class of normal SNe la with the relatively 
large amount of ^®Ni produced by the explosion of cold 
Chandrasekhar-mass white dwarfs. The most promising 
theoretical framework for over two decades has been the 
deflagration-to-detonation transition (DDT), in which a 
turbulent flame ignited within the convective core of a 
Chandrasekhar-mass white dwarf leads to slow, subsonic 
burning, and a pre-expansion of the white dwarf itself 
(Khokhlov 1991; Gamezo et al. 2005). The DDT conjec¬ 
tures that the flame undergoes a subsequent transition 
to a detonation, with the material in the pre-expanding 
star producing a normal la event. However, as we explain 
here, recent theoretical work on the convective simmer¬ 
ing phase within Chandrasekhar-mass white dwarfs is in¬ 
consistent with the large amount of burning required by 
the standard DDT model for normal SNe la during the 
deflagration phase in the majority of ignitions. 

In this paper, we present the case that single¬ 
degenerate SNe la generally lack a vigorous deflagration 
phase, and preferentially result in overluminous SNe la. 
This hypothesis is motivated by the need to reconcile 
numerous recent observational results, ranging from the 
DTD (Totani et al. 2008; Graur et al. 2011; Maoz & 
Badenes 2010; Maoz et al. 2010, 2012; Graur et al. 2014), 
the absence of evidence for companions (Li et al. 2011; 
Bloom et al. 2012) and ex-companions (Maoz & Man- 
nucci 2008; Gonzalez Hernandez et al. 2012; Kerzendorf 
et al. 2012; Schaefer & Pagnotta 2012; Edwards et al. 
2012; Kerzendorf et al. 2014), the absence of H in the 
nebular phase (Leonard 2007), and the absence of radio 
and X-ray emission (Li et al. 2011; Brown et al. 2012; 
Horesh et al. 2012; Bloom et al. 2012; Margutti et al. 
2012; Chomiuk et al. 2012) around peak brightness, with 


recent theoretical advances. In particular, recent theo¬ 
retical work has, for the first time, clarified the inher¬ 
ently stochastic nature of carbon ignition in the convec¬ 
tive cores of Ghandrasekhar-mass white dwarfs. The lit¬ 
erature of explosion models has focused primarily either 
upon precisely-centered and significantly off-centered ig¬ 
nitions, with widely ranging outcomes dependent upon 
the ignition - from subluminous 2002cx-like events (Jor¬ 
dan et al. 2012b; Kromer et al. 2013), though overlumi¬ 
nous SNe la events (Plewa et al. 2004; Jordan et al. 2008; 
Meakin et al. 2009). Here, we develop an analytic crite¬ 
rion that separates between buoyancy-driven and central 
ignitions. Criteria for central ignition have been explored 
previously in the literature (Zingale & Dursi 2007; As- 
pden et al. 2011; Malone et al. 2014) using numerical 
solutions to semi-analytic models and fully multidimen¬ 
sional simulations. Here, we provide a new, elementary 
derivation of the critical length scale within which cen¬ 
tral ignitions may arise, which highlights the essential 
competition between the physical processes of burning 
and buoyancy. We further demonstrate that the major¬ 
ity of ignition events within the single-degenerate chan¬ 
nel lead to buoyancy-driven turbulent deflagration, with 
little deflagration energy release and pre-expansion, and 
argue these will be followed by detonation and conse¬ 
quently an overluminous SNe la. We place bounds on the 
number of normal and subluminous SNe la events orig¬ 
inating from the single-degenerate channel, and demon¬ 
strate that these can only account for a small fraction 
of all single-degenerate events. We argue that normal 
and subluminous events must primarily arise through the 
double-degenerate and helium-donor channels. 

These findings clarify a number of long-standing prob¬ 
lems, including the rate problem of the single-degenerate 
channel, the absence of observed stellar progenitors and 
ex-companions, and the absence of hydrogen in the neb¬ 
ular phase of normal supernovae. They also connect 
the particular class of overrluminous SNe la to stellar 
progenitors, including recurrent novae and supersoft X- 
ray sources, and strengthen the association of single¬ 
degenerate channel SNe la to the subclass of SN la-CSM 
events like PTFllkx (Dilday et al. 2012). Additionally, 
they shed additional light on why overluminous super¬ 
novae primarily originate in late-type galaxies. Our re¬ 
sults offer guidance to observers to look towards overlu¬ 
minous events for the most likely occurrence of evidence 
either to support or constrain the single-degenerate chan¬ 
nel. 

In §2, we review the theory of single-degenerate 
Chandrasekhar-mass white dwarfs, from ignition through 
detonation. In §3, we delve into the ignition problem in 
depth, and develop an analytic criterion that separates 
between buoyancy-driven and central ignitions. In §4, we 
compare the rates both of supersoft X-ray sources (SSSs), 
as well as binary population synthesis models, against 
those of overluminous la events. In §5, we conclude with 
the observational implications of this work, and suggest 
specific observations which can either directly support 
or constrain the single-degenerate channel as a primary 
contributor to overluminous SNe la events. 

2. SINGLE-DEGENERATE CHANNEL 

The single-degenerate channel has been the most 
thoroughly-explored SNe la model framework over the 
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last several decades (Whelan & Iben 1973). It begins 
with an intermediate mass primary star between 3 Mq 
- 8 Mq and a companion. The primary evolves more 
rapidly, and transfers mass to the companion, leading 
to a common envelope phase of evolution. The primary 
forms a 0.6 - 1.2 Mq C/0 WD, which can then reach 
the Chandrasekhar mass by hydrogen-rich mass trans¬ 
fer from the non-degenerate companion. Modeling this 
sequence of events using theoretical binary population 
synthesis (BPS) codes demonstrates the broad feasibility 
of the channel, although the predicted rates are subject 
to considerable uncertainties, which we will review in §4. 
Observationally, the accreting white dwarfs are visible as 
recurrent novae and supersoft X-ray sources. The more 
rapidly-accreting supersoft sources accumulate mass up 
to the Chandrasekhar mass limit, and contrary to theo¬ 
retical expectations, recent observational evidence of the 
mass functions of pre-CV and post-CV WDs seems to 
suggest that recurrent novae may also accumulate mass 
(Zorotovic et al. 2011). 

The central conditions of accreting white dwarfs reach 
carbon ignition as the white dwarf approaches the Chan¬ 
drasekhar mass limit (Nomoto et al. 1984). A vari¬ 
ety of explosion models have been considered, rang¬ 
ing from pure deflagrations (Ropke et al. 2007a; Ma 
et al. 2013), through deflagration-to-detonation transi¬ 
tions (DDTs) (Gamezo et al. 2005; Ropke & Niemeyer 
2007), and gravitationally-confined detonations (CCDs) 
(Plewa et al. 2004; Townsley et al. 2007; Jordan et al. 
2008; Meakin et al. 2009; Jordan et al. 2012a,b). 

It has slowly become clear that, depending on the 
nature of ignition, the single-degenerate channel pre¬ 
dicts a wide range of deflagration energy yields. The 
greater the deflagration energy release, the greater the 
pre-expansion, and the lesser the yield of ®®Ni during the 
following detonation phase. Thus, the single-degenerate 
channel predicts a range of outcomes, from overlumi- 
nous events like SN 1991T (Jordan et al. 2008), through 
normal la events (Gamezo et al. 2005; Jordan et al. 
2012a), and even, with highly vigorous deflagration burn¬ 
ing and in the absence of any detonation, subluminous 
SN 2002cx-like events (Jordan et al. 2012b; Kromer et al. 
2013). The energetic yield during the deflagration phase 
is directly tied to the ignition problem, which until very 
recently was considered a largely unsolved problem, with 
theorists addressing the uncertainty by widely varying 
the number and offset location of ignition points in mul¬ 
tidimensional models of the deflagration and detonation 
phases. However, a monumental computational cam¬ 
paign, primarily by Zingale, Woosley, Bell and collab¬ 
orators (Zingale et al. 2011; Nonaka et al. 2012; Malone 
et al. 2014) has revealed that the simmering phase typ¬ 
ically yields a single ignition point, with a distribution 
of offset ignition locations, and an expectation value of 
^ 50 km. 

The distribution of offset locations reflects the inher¬ 
ently stochastic nature of turbulent convection. This ir¬ 
reducible stochasticity of thermonuclear ignition directly 
results in an inherently probabilistic range of outcomes 
for the single degenerate channel. Sufficiently large offset 
radii lead to buoyancy-driven ignition and deflagration, 
but precisely how offset must an ignition point be to 
lead to this outcome? This is the key question which we 
ask and address in the next section, where we quantify 


the probability of a central ignition by deriving a charac¬ 
teristic length scale separating buoyancy-driven ignitions 
from central ignition events. 

3. IGNITION AND DETONATION 

3.1. Central and Offset Ignition 

A single buoyancy-driven ignition point burns a small 
fraction of the white dwarf during the deflagration phase 
(Plewa et al. 2004; Townsley et al. 2007; Ropke et al. 
2007b; Jordan et al. 2008; Meakin et al. 2009; Malone 
et al. 2014), and consequently leads to small amounts 
of pre-expansion, consistent with an overluminous SNe 
la event. In contrast, central ignition points linger 
near the central region of the star, leading to greater 
pre-expansion, typically producing ~ 0.5Mq of ®®Ni 
burnt during the deflagration phase (Nomoto et al. 1984; 
Gamezo et al. 2005; Ma et al. 2013). This greater amount 
of burning during the deflagration phase may yield a sub- 
luminous SNe la by itself (Jordan et al. 2012a; Kromer 
et al. 2013), or a normal SNe if followed by a detonation, 
in either the context of a DDT or CCD (Gamezo et al. 
2005; Meakin et al. 2009; Jordan et al. 2012b). Hence, for 
single-bubble ignitions, a key question is which ignition 
radii lead to buoyany-driven ignitions, and which lead to 
central ignitions. The answer to this question is crucially 
connected to the amount of pre-expansion which the WD 
experiences during the deflagration phase, which in turn 
determines whether a DDT or GCD is possible, as we 
will see in section 3.2. 

We now develop an analytic criterion which demarcates 
the boundary between buoyancy-driven and central igni¬ 
tions. We assume a single ignition bubble, with offset ra¬ 
dius To from the center of the white dwarf, and radius R 
- see figure 1. We assume that the ignition size is small in 
comparison to the offset radius - R Tq. A crucial criti¬ 
cal length scale in the development of the flame bubble is 
the fire-polishing scale Afp = AnSf /Ag. Perturbations to 
the flame bubble on length scales below Afp are polished 
away by the action of the flame, whereas larger-scale per¬ 
turbations are subject to the Rayleigh-Taylor instability 
(Timmes & Woosley 1992). For typical ignition condi¬ 
tions, the bubble radius R is much less than the fire¬ 
polishing scale Afp and as a result, is necessarily spheri¬ 
cal upon ignition. Furthermore, the bubble radius is also 
much smaller than the integral length scale L ~ 100 km 
of the convective turbulent background. Consequently, 
the turbulence associated with convection is also negligi¬ 
ble on the scale of the bubble during its initial growth. 

While the neglect of both the non-spherical distortion 
of the bubble and of turbulence is an excellent approx¬ 
imation initially, they begin to break down as the bub¬ 
ble burns outward (Malone et al. 2014). However, for 
the purpose of determining which ignition points are 
buoyancy-driven, we will be interested in only the ini¬ 
tial growth of the bubble, during which both approxima¬ 
tions continue to hold valid throughout. Furthermore, 
numerical simulations demonstrate that the bubbles are 
typically ignited within a convective plume with a net ra¬ 
dial velocity, though the magnitude of the typical plume 
velocity (^ 20 km/s) is relatively small in comparison to 
the laminar flame speed Si ~ 100 km/s, and may also be 
neglected to a Hrst approximation. Further, we demon¬ 
strate below that the neglect of the plume velocity places 
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Figure 1. A schematic diagram depicting the geometry of the 
critical case of a bubble, ignited at a distance rg from the center 
of the star, which just barely begins to burn through the center of 
the star in the time required for it to buoyantly rise an additional 
vertical distance r'. As shown, the radius H of the critical case 
bubble at the instant it begins to burn through the center of the 
star is H = rg + r'. 


a conservative upper-bound on the possible central igni¬ 
tion events. Therefore, we will consider the initial growth 
of the bubble on a purely stationary background. 

The criterion for central ignition is simply expressed; 
the bubble radius i?(t) must exceed the offset radius r(t) 
at some time t - that is, R{t) > r{t). A central ignition 
burns through the center of the star prior to buoyantly 
rising an additional vertical distance r'. In the appendix 
material, we consider the combination of buoyancy, drag 
forces, and the added mass effect on a self-similar spher¬ 
ical flame bubble that develops from a point-like ignition 
at t = 0 - that is, R{t = 0) = 0. There, we demonstrate 
that the combination of these effects result in an effec¬ 
tive gravitational acceleration geg. That is, the leading 
terms of a power-series solution of the equation of motion 
of the bubble, an ordinary differential equation describ¬ 
ing the radial location r, equation (A8), are given by 
r(t) = ro -b (( 7 eff/ 2 )t^ at early times after ignition, with 
geB is given by 

^ (3-A) 

Here A is the Atwood number, a dimensionless quantity 
describing the density contrast across the bubble inter¬ 
face : 

A = PlAlfA (2) 

Pa+ Pf 

where Pf is the density of the fuel, and pa is the density 
of ash, and the effective gravitational acceleration g^s is 
taken at the initial radius rg. The density of fuel and 
ash are similar under carbon burning at central densi¬ 
ties typical of Chandrasekhar-mass white dwarfs, lead¬ 
ing to small Atwood number values A ~ 0.09 (Timmes 
& Woosley 1992). As shown in the appendix, in equation 
(A5), in the absence of both drag and the added mass 
effect, the effective buoyant acceleration would simply be 

Ap , , 

geS= 9 — (3) 

P 

where Ap/p = 2AI{1 + A). Consequently, the effective 
gravity for a low Atwood number bubble (A <C 1) at early 
times, including the effects of drag and added mass, is 
very nearly 1 /6 of the buoyancy-driven value that would 


be obtained in the absence of both drag and the added- 
mass effect. 

The bubble achieves a buoyant velocity v = y/2gesR 
after rising the radial distance r'. In the critical case, the 
buoyant velocity v will exactly equal the radial expansion 
velocity of the bubble [(1 -I- A)/(l — A)]5'; (equation Al), 
and the bubble just barely burns through the center as 
it reaches the offset radius rg -|- r'. Hence, 


! 

r 


i+A y s‘f 

1 — A / 2peff 


(4) 


Furthermore, the critical condition that the radius R 
of the bubb le reach es the center of the star in this same 
time, t = ^2r'/Peff yields 


i? = rg -I- r' 


1 + A y 

1 — A / peff 


( 5 ) 


Consequently, in the critical case, the additional distance 
the bubble must rise to burn through the center of the 
star is precisely equal to its offset radius: r' = rg. There¬ 
fore, a buoyancy-driven ignition must exceed a critical 
offset radius rcrit 


ro > rcrit = 


/1 + AA 

1 ^ 1 

n+A\ 

U-aJ 

' 2peff StT ' 

U-aJ 


Afp 


( 6 ) 


This result demonstrates that the critical offset radius for 
a buoyancy-driven ignition event is proportional to, but 
signihcantly smaller than the hre-polishing scale, with 
the laminar flame speed and effective gravity taken at 
the ignition location. 

The competition between laminar burning and buoy¬ 
ancy is intrinsically linked through the fire-polishing 
scale Afp. In particular, because the fire-polishing scale 
is inversely proportional to the effective gravitational ac¬ 
celeration, near the center of the star, where g^s 0, 
the laminar burning timescale is shorter than the buoy¬ 
ancy timescale of the bubble, and bubbles ignited there 
must burn through the center of the star. Conversely, at 
larger offset radii, as the effective gravitational accelera¬ 
tion increases, the fire-polishing scale Afp drops, and the 
buoyancy timescale is shorter than the laminar burning 
timescale. The critical radius rcrit defined by equation 6 
establishes the boundary between these two regimes. 

The local effective gravitational acceleration may be 
expressed as peS = (dTr/SjGeffPcrg, with G^b = [A/(3 — 
A)]G the effective gravitational coupling for the bubble: 


^crit — 


1-f A 
1 - A 


STrGeffP, 


-Si 


(7) 


Further insight into the physics of the buoyancy-driven 
regime can be gathered by noting that buoyancy-driven 
ignition occurs when the initial offset radius rg exceeds a 
value which is proportional to the product of the laminar 
flame speed and the central buoyant dynamical time : 


^crit — 

TT 


l-f A 
1-A 

= 19 km 


Sltgy 


Si 

100 km/s 


^.09y/^ ^ 2.2 X 10^ gm/cm^ y/^ 
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Here the dynamical time is defined as the buoyancy dy¬ 
namical timescale tdyn = a/Stt/( 32060/3), based upon the 
effective gravitational acceleration, and the the fiducial 
scales are determined for a 50/50 C/0 WD (Timmes & 
Woosley 1992). 

The critical radius rcrit is very close to the center 
of the WD, and highlights the importance of an ana¬ 
lytic derivation. In figure 2, we plot the analytic solu¬ 
tion trajectory of bubble radius versus offset radius for 
a buoyancy-driven ignition, a central ignition, and the 
critical ignition case. For comparison, we also overplot 
a full numerical integration of equations (Al) and (A3) 
on a fully-stratified white dwarf background, for which 
we determined a critical ignition radius rcrit = 24.1 km, 
and show very good agreement with the critical solu¬ 
tion trajectory. Large-scale three-dimensional full-star 
numerical simulations typically have spatial resolutions 
ranging from 4-8 km, and therefore only marginally re¬ 
solve the critical transition radius at best. Some compar¬ 
isons, can however, be made between our analytic result 
and those obtained from previous analytic modeling and 
three-dimensional numerical simulations. Aspden et al. 
(2011) developed a semi-analytic model very similar to 
the analytic one presented in the appendix, though in¬ 
cluding both an entrainment term for the bubble growth 
as well as two free parameters which were then fit to 
data from three-dimensional models of the ignition and 
early evolution of the flame bubble, and then solved nu¬ 
merically on a full-star background. Aspden et al. (2011) 
present evidence for a transition to central ignition some¬ 
what interior to Tq = 16 km, though it should be noted 
that they used a bubble with an initial radius of 10.2 km, 
which is comparable to the critical radius, and adopted a 
somewhat higher central density 2.55 x 10® g cm“^, than 
assumed here. Zingale & Dursi (2007) numerically find 
?’crit ~ 23.5 km, using small initial bubble radii like our 
self-similar solution, and in very good agreement with our 
own numerical solution. However, they assume a higher 
central density of 2.6 x 10® g cm“® and employ a fit for the 
laminar flame speed that is modified to improve its ac¬ 
curacy at low densities compared to the fit from Timmes 
& Woosley (1992) that we use, in addition to assuming a 
lower carbon mass fraction of 0.3. A more recent study 
by Malone et al. (2014) finds a critical radius of ^ 10 
km over a range of point-like initial bubbles. Similarly, a 
variety of other two and three-dimensional full-star sim¬ 
ulations find a significant and sharp enhancement in the 
burning during the deflagration phase for rg ^ 20 km — 
see e.g. Plewa (2007) and Meakin et al. (2009), though 
the resolutions in these models is relatively coarse, and 
begin with bubbles with radial extents comparable to the 
critical radius. Overall, our analytic result is, however, 
consistent with these previous findings. 

Numerical simulations of the convective phase of 
Chandrasekhar-mass white dwarfs have revealed the dis¬ 
tribution of hot spots within the core. Zingale et al. 
(2011) and Nonaka et al. (2012) have examined the dis¬ 
tribution of hot spots leading up to ignition. Due the 
strongly nonlinear character of the ignition problem, only 
one of these hot spots ignites and runs away. We note 
that, due to the significant computational cost of these 
three-dimensional numerical simulations of convection 
and ignition, only a single progenitor model has been ex- 



Figure 2. A plot of bubble radius R versus offset radius ro. The 
dashed line R = r is the critical case separating central ignitions 
from buoyancy-driven ignitions. The hatched region of the figure 
marks central ignitions, with R > r; buoyancy-driven ignitions 
lie in the portion of the figure with r < R. The vertical dashed 
line shows the critical offset radius, and the solid lines denote the 
trajectories of the analytic solution described in the text, with 
the thick solid line showing the critical case of a bubble which 
barely burns through the center of the star. For comparison, a full 
numerical integration of a critical bubble in a fully-stratified WD, 
with an initial critical offset = 24.1 km is shown in the dotted 
curve. 

amined in detail to date, and additional models will be 
needed to further elucidate the nature of the ignitions. 
However, we consider this distribution of hot spots as 
representative of the temperature distribution of all pos¬ 
sible ignition points within the white dwarf, and therefore 
representative of the statistical distribution which would 
be obtained for the final ignition points over an ensemble 
of simulations. 

Conservatively, from direct inspection of their figure 7, 
we find that no greater than 2 ± 1% of the ignition hot 
spots considered by Nonaka et al. (2012) burn through 
the center of the star, and lead to central ignitions, with 
a commensurately-larger energy release during the de¬ 
flagration phase. The remaining hot spots will lead to 
buoyancy-driven ignition events. We emphasize that 
these estimates are conservative upper-bounds on the 
likelihood of a central ignition, due to our neglect of the 
generally upward plume velocities which the hot spots are 
ignited within, and to the enhanced effect of buoyancy in 
the numerical simulations over analytic models (Malone 
et al. 2014). Thus, the overwhelmingly likely outcome 
of ignition in a Chandrasekhar-mass white dwarf is a 
buoyancy-driven event. As we demonstrate in the next 
section, these buoyancy-driven ignition events are in turn 
generally expected to lead to overluminous SNe la. 

3.2. Detonation 

Pure deflagration models (Ropke et al. 2007a; Ma 
et al. 2013) have numerous deficiencies when compared 
against observations, including significant amounts of un¬ 
burned carbon and oxygen, which can be rectified if the 
white dwarf undergoes a subsequent detonation. Hence, 
for many years, the leading single-degenerate channel 
model has been the deflagration-to-detonation transition 
(DDT) (Khokhlov 1991; Gamezo et al. 2005; Ropke & 
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Niemeyer 2007). However, the DDT model crucially rests 
on a significant amount of nuclear burning during the de¬ 
flagration phase, prior to the subsequent detonation, in 
order to pre-expand the WD and achieve a nucleosyn- 
thetic yield of ®®Ni typical of normal brightness SNe la. 
A key consequence of the central ignition criterion, equa¬ 
tion (6), is that such a vigorous deflagration phase is 
highly unlikely in the majority of single-degenerate SNe 
la ignition events. Instead, the overwhelmingly-likely 
outcome of ignition is a buoyancy-driven flame bubble 
with a relatively small deflagration energy yield. 

The subsequent evolution of a buoyancy-driven flame 
bubble has been explored in numerous calculations in 
the literature. In a gravitationally-confined detonation 
(GCD), the small amount of nuclear energy released dur¬ 
ing the deflagration phase leads to a bound white dwarf, 
and the hot bubble breaks out of the bound white dwarf 
surface. In a GCD, the white dwarf undergoes a det¬ 
onation as the hot ash remaining from the bubble is 
ram-driven into high densities (Plewa et al. 2004; Towns- 
ley et al. 2007; Jordan et al. 2008). We note that these 
findings were not initially universally agreed-upon, with 
another group using an enhanced subgrid model for tur¬ 
bulent combustion found a much greater degree of pre¬ 
expansion, leading to a detonation in 2D but not 3D 
(Ropke et al. 2007b). However, more recent work, both 
semi-analytic and numerical, has demonstrated the de¬ 
flagration energy release is established by large-scale hy¬ 
drodynamic entrainment of the bubble, and is less sensi¬ 
tive to the turbulent flame model adopted (Aspden et al. 
2011). Additionally, a third group (Malone et al. 2014) 
modeling the convective phase though ignition and bub¬ 
ble breakout has obtained deflagration results very sim¬ 
ilar to that of the classic GCD model of Jordan et al. 
(2008), using a completely independent code. 

There are several caveats to consider when evaluating 
this theoretical work. Firstly, our physical and computa¬ 
tional description of turbulent deflagration is still limited. 
Much of the burning during the deflagration phase occurs 
on length scales unresolved in multi-dimensional numer¬ 
ical simulations. Several approaches have been advanced 
to address with this issue. While there are some differ¬ 
ences in the amount of burning in the deflagration pre¬ 
dicted by these approaches, they generally agree in that 
the amount of nuclear energy released for a single bubble 
ignition is small in comparison to the gravitational bind¬ 
ing energy of the WD. Only a significant revision of our 
basic understanding of turbulent combustion would fun¬ 
damentally alter this outcome. Moreover, recent analytic 
and semi-analytic work on burning bubbles has clarified 
that the large-scale hydrodynamical entrainment plays 
a more important role than the choice of the turbulent 
flame speed (Aspden et al. 2011). 

Secondly, three-dimensional simulations of both the 
convective phase and the subsequent deflagration and 
detonation phases typically rely upon a standard G/0 
WD model, with 50/50 G/0 composition, and masses 
1.36 - 1.38 Mq. However, the deflagration phase is sensi¬ 
tive to the value chosen for this mass, primarily through 
the nonlinear dependence of the laminar flame speed to 
the density, leading to commensurately greater flame pol¬ 
ishing at the higher central densities for WD progeni¬ 
tors approaching the Ghandrasekhar mass (Krueger et al. 
2012; Lamb et al. 2014). Higher laminar flame speeds 


polish the flame bubble surface, and consequently pro¬ 
duce less burning than a lower laminar flame speed. Even 
a very small relative change in the WD mass, of 1% or 
less, at ignition can produce a significant change in the 
energy yield in the deflagration phase. Thus it is possible 
to produce a range of outcomes, even including normal 
SNe la, by fine-tuning of the chosen progenitor WD mass 
very close to the Ghandrasekhar mass. We view such a 
fine-tuning as unlikely, since it would tend to produce 
a relatively large rate of normal single-degenerate SNe 
la, in tension with observations. However, this is an im¬ 
portant issue which needs to be understood better, since 
it may represent the dominant theoretical uncertainty of 
theoretical modeling of the single-degenerate channel. 

Additionally, a buoyancy-driven flame bubble may un¬ 
dergo a transition to detonation before it even breaks out 
from the white dwarf surface. Such a transition to det¬ 
onation would differ significantly from the classic DDT 
scenario, however, due to the lower deflagration energy 
release, with the subsequent outcome still inevitably an 
overluminous event. Lastly, the first bubble may not pro¬ 
duce an effective explosion, and may instead be followed 
by a series of successive bubble ignitions, only one of 
which is required to detonate the white dwarf, as specu¬ 
lated recently by Malone et al. (2014). 

In summary, the bulk of the numerical evidence demon¬ 
strates that a detonation must follow buoyancy-driven 
events. The GGD scenario is the most thoroughly- 
explored model, and the most likely result of a buoyancy- 
driven ignition. However, any of the variants described 
above will generally produce overluminous events, pro¬ 
vided that a detonation follows the ignition of one or a 
small number of ignition events. 

4. RATES, HOST GALAXY ENVIRONMENTS, AND 
DELAY-TIME DISTRIBUTIONS 

A long-standing problem in all major supernova chan¬ 
nels is how to reconcile the observed SNe la rates with 
those predicted both from observations of likely stellar 
progenitors, as well as theoretical binary population syn¬ 
thesis models. As mentioned previously, estimates of the 
rates of overluminous SN 1991T-like events range from 
4% - 20% of the overall SNe la rate (Foley et al. 2009; 
Li et al. 2001). In this section, we demonstrate that the 
need to account for only overluminous SNe la, as opposed 
to the entire population of SNe la, significantly alleviates 
much of the tension on single-degenerate rates. We also 
discuss the delay-time distributions of single-degenerate 
models, and compare these to the observed galactic host 
environments of overluminous SNe la. 

Supersoft X-ray sources have long been thought to be 
possible SNe la progenitors - e.g. van den Heuvel et al. 
(1992). However, the number of known supersoft X- 
ray sources (Di Stefano 2010) and the X-ray background 
(Gilfanov & Bogdan 2010) are discrepant with observed 
SNe la rates by 1 - 2 orders of magnitude. 

Rates derived from theoretical binary population syn¬ 
thesis models span a very wide range, depending on 
model assumptions. In particular, the derived rates are 
highly sensitive to the prescription of mass retention 
adopted, and to a lesser extent also the assumptions re¬ 
garding common envelope evolution and mass transfer 
(Hours et al. 2013; Glaeys et al. 2014). Recent work 
suggests that the wide range of predictions made by dif- 
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ferent BPS models employed by different groups can be 
connected back to varying model assumptions (Toonen 
et al. 2014). 

Bearing in mind these very significant uncertainties, 
we can compare the SNe la rates determined from BPS 
models with those observed both for normal and over- 
luminous SNe la. The time-integrated BPS rates in the 
single-degenerate channel range from .35 x 10 “"^Mq^ — 
1.3 X 10“^Mq^ (Bours et al. 2013). In contrast, the ob¬ 
served time-integrated SNe la rates are in the range of 
4 X 10-^Mq ^ - 26 X IO-^Mq ^ (Maoz et al. 2012; Maoz 
& Mannucci 2012; Perrett et al. 2012; Graur & Maoz 
2013). In other words, current BPS models suggest that 
the single-degenerate channel accounts for only 1% - 30% 
of the total SNe la rate. However, this range spans the 
range of estimates of overluminous SNe la. 

In contrast to buoyant ignitions, central ignitions in 
single-degenerate SNe la lead to significant amounts of 
deflagration energy release, and could possibly produce a 
range of outcomes, from failed 2002cx-like SNe le events 
through normal events. These are, however, less likely 
outcomes given the significantly smaller probability of 
ignition within the central ignition region. We now es¬ 
timate the single degenerate contribution to normal and 
failed SNe la, with two key assumptions. First, we posit 
that the fraction of normal and subluminous SNe la 
events stemming from the single-degenerate channel is 
equal to the fraction of central ignitions compared to the 
total. Secondly, we further assume that the overall pop¬ 
ulation of overluminous events primarily originates from 
the single-degenerate channel. That is, if other channels, 
including the double-degenerate channel, do contribute 
to the population of overluminous SNe la, as we discuss 
below in §5, we assume that their contribution is negli¬ 
gible compared to that of single-degenerates. We place a 
conservative upper bound on the single-degenerate con¬ 
tribution by taking the fraction of central ignitions to be 
5% (including a 3 cr variance) and the observed fraction 
of overluminous SNe la to be 20%. Consequently, we 
infer the contribution of the single degenerate channel 
to both normal SNe la and failed 2002cx events to be 
no greater than 1% of the overall SNe la rate. The low 
probability of central ignition generally implies the con¬ 
tribution of single-degenerates to normal and failed SNe 
la events must be small. In particular, if other channels 
also contribute signihcantly to the overluminous SNe la 
rate, the single-degenerate contribution to normal and 
failed SNe la events will be even smaller than our con¬ 
servative upper-bound. We thus suggest that the progen¬ 
itors of normal and failed 2002cx events must primarily 
originate from other channels - most likely the double¬ 
degenerate and helium donor channels. 

Many studies have confirmed that overluminous SNe 
la events arise primarily in late-type galaxies Howell 
(2001), and that late-type galaxy SNe la are overlumi¬ 
nous in comparison to those which arise early-type galax¬ 
ies (Hamuy et al. 2000; Sullivan et al. 2006; Lampeitl 
et al. 2010; Gupta et al. 2011; Johansson et al. 2013; 
Childress et al. 2013). The late-type galaxy host en¬ 
vironment of overluminous SNe la can be explained in 
the context of the single-degenerate channel by noting 
that the single-degenerate DTD is relatively prompt in 
comparison to the double-degenerate channel, with de¬ 


lay times starting at 40 - 100 Myr, depending on model 
assumptions (Nomoto et al. 2007; Ruiter et al. 2009; Too¬ 
nen et al. 2014). Moreover, the single-degenerate channel 
DTD is distinguished from the double-degenerate chan¬ 
nel in the need to have a main-sequence or red giant 
donor; consequently the single-degenerate DTD exhibits 
a sharp cutoff at 2 - 3 Gyr, independent of model as¬ 
sumptions. Hence, single-degenerate SNe la are natu¬ 
rally expected to be present in relatively young stellar 
populations, and absent from very old (> 2 - 3 Gyr) 
stellar populations, as has long been noted - see e.g. 
della Valle & Livio (1994). However, we also note that 
massive C/0 primaries are naturally more luminous in 
both the double-degenerate and double-detonation mod¬ 
els as well. Such massive C/0 WD primaries arise from 
more massive intermediate-mass stellar progenitors, with 
prompt delay times (< 100 Myr) (Ruiter et al. 2013), 
and therefore also naturally arise in late-type galaxies. 
Hence the correlation between peak brightness and host 
galaxy environment is one which is, broadly speaking, 
expected for single-degenerate SNe la as well as double¬ 
degenerate and double-detonation SNe la. Consequently, 
in order to better understand the correlation between the 
host environments and luminosities of SNe la, additional 
observational constraints which directly test the single¬ 
degenerate channel are needed, as we discuss in the dis¬ 
cussion in §5. 

In summary, if the single-degenerate channel must ac¬ 
count only for overluminous SNe la, and not the popu¬ 
lation of normal SNe la as a whole, the estimated rates 
derived both from the known population of supersoft X- 
ray sources and from theoretical BPS models are much 
more closely aligned. While much work remains to be 
done, particularly with regards to improving the mass 
retention efficiency physics of the BPS models, the over¬ 
all consistency between the rates of the single-degenerate 
sources and overluminous SNe la rates motivates refined 
observations to test this connection further. 

5. DISCUSSION 

Sim et al. (2013) have demonstrated that synthetic 
spectra of single and few-bubble offset ignition simu¬ 
lations are nominally good matches to SN 1991T-like 
SNe la events, according to the standard criterion of 
a sufficiently high spectral cross-correlation coefficient 
(Blondin & Tonry 2007). The synthetic spectra support 
the view that single-degenerate supernovae are promis¬ 
ing candidates for SN 1991T-like events. Because the 
deflagration stage is most particularly sensitive to the 
central density of the progenitor, further exploring this 
model sensitivity may help clarify the quality of the 
cross-correlation of models with observations, and pro¬ 
vide theoretical spectral templates to compare against 
observations. 

In principle, the presence of a main-sequence or red- 
giant companion star is detectable in early-time light 
curve distortions arising from the shocked gas of the com¬ 
panion under favorable viewing angles (Kasen 2010). Ex¬ 
tensive light curve data from the Sloan Digital Sky Sur¬ 
vey and the Supernova Legacy Survey (Hayden et al. 
2010; Bianco et al. 2011) failed to produce evidence for 
such distortions, and constrained the fraction of red giant 
companions to be 10% - 20% of the overall population 
of SNe la - within the range expected for overluminous 


SNe la. Further exploration of these early light curve sig¬ 
natures focusing upon overluminous SNe la events may 
be an interesting avenue for future work. However, the 
recent identification of the diversity of the distribution 
of 56Ni in the outermost ejecta of SNe la, reflected in 
a range of early light curve profiles (Firth et al. 2015) 
will make it more challenging to characterize the possi¬ 
ble presence of shocked gas from the companions in the 
data. 

The presence of any significant amount of hydrogen 
in the nebular phase of a SNe la spectrum would be a 
“smoking gun” piece of evidence of stripped hydrogen 
from the companion star in favor of the single-degenerate 
channel, but has so far eluded detection. Tight obser¬ 
vational constraints have been placed on the hydrogen 
abundance in SNe, as determined from the absence of 
any observable Ha emission lines in the nebular phase 
spectra of normal la. Leonard (2007) constrained the 
amount of H in the normal SN la 2005cf and the slightly 
subluminous SN la 2005am to be less than 
Shappee et al. (2013b), Lundqvist et al. (2015), and Gra¬ 
ham et al. (2015) applied similar methods to place the 
tightest-known constraints on the amount of H in SN 
201 Ife, at < 3 X IO^^Mq. Although theoretical caveats 
remain to be investigated that may serve to obscure Ha, 
even at late times - see e.g., Leonard (2007) - current 
observational limits on SN 201 Ife and other normally- 
bright SNe la are well beneath the range predicted for 
nearly all types of red giant and main sequence, including 
helium main sequence companions (Marietta et al. 2000; 
Pakmor et al. 2008; Pan et al. 2013, 2014), with the pos¬ 
sible exception of M-dwarfs (Wheeler 2012). Crucially, 
however, we note that the most important constraints on 
Ha in the nebular phase of SNe la are all derived from 
normal SNe la. The presence of Ha in overluminous SNe 
la remain unconstrained by data of comparable quality. 

Consequently, one of the most important tests of the 
contribution of single-degenerate supernovae to overlu¬ 
minous SNe la will be a deep probe for Ha in the nebu¬ 
lar optical spectra of low-redshift overluminous SNe la. 
Additional recent work suggests Pa in the near-infrared 
post-maximum (Maeda et al. 2014) may also be a possi¬ 
ble strong discriminant of the companion. Ongoing sur¬ 
veys including ASASSN will identify a number of low- 
redshift SN 1991T-like SNe la which will enable a com¬ 
prehensive survey in a number of events. Already sev¬ 
eral good candidates, including ASASSN-141w, PESSTO 
ESO 154-GlO, and ASASSN-15bc at z < .04 have re¬ 
cently been discovered, and will be prime candidates for 
nebular-phase follow-up in the near future. 

The relatively rare class of SNe la-CSM demonstrates 
evidence for Ha emission. SNe la-CSM originate exclu¬ 
sively in late-type galaxies and irregulars and are overlu¬ 
minous, with —21.3 < Mfi < —19 peak mag (Silverman 
et al. 2013). Moreover, recent evidence points to an as¬ 
sociation of SNe la-CSM and overluminous SNe 1991T- 
like events (Leloudas et al. 2015). In contrast, multi¬ 
wavelength observations, including upper-bounds of ra¬ 
dio and X-ray emission from the normal SN la 201 Ife, 
place tight constraints on the possible stellar compan¬ 
ions and the circumstellar environment around normal 
SNe la (Li et al. 2011; Brown et al. 2012; Horesh et al. 
2012; Bloom et al. 2012; Margutti et al. 2012; Chomiuk 
et al. 2012). Similar observations, carried out for an over- 


luminous SNe la event at a comparable distance, are the 
most promising site of the first multi-wavelength signa¬ 
ture from a SN la. 

Another strong piece of evidence in favor of the single¬ 
degenerate channel would be the detection of the com¬ 
panion star in a Type la supernova remnant. How¬ 
ever, until very recently, there has been a general lack of 
supporting evidence of these so-called “ex-companion” 
stars. In particular, observers have found an absence 
of ex-companions in SNRs SN2006dd and SN2006mr in 
NGC 1316 (Maoz & Mannucci 2008), SN 1006 (Gonzalez 
Hernandez et al. 2012; Kerzendorf et al. 2012), SNR 
0509-67.5 (Schaefer & Pagnotta 2012), SNR 0519-69.0 
(Edwards et al. 2012), and SN 1604 Kepler (Kerzendorf 
et al. 2014). A candidate ex-companion in SN 1572 Ty¬ 
cho has been identified (Ruiz-Lapuente et al. 2004), but 
has been questioned as our knowledge of the properties of 
ex-companions improved (Shappee et al. 2013a; Pan et al. 
2014). Very recent work has revealed a number of possi¬ 
ble ex-companions in the SNe la remnants SNR 0505- 
67.9 and SNR 0509-68.7 (Pagnotta & Schaefer 2015). 
While the sample size of investigated remnants is so far 
very small, if one or both of the ex-companions in the 
remnants 0505-67.9 and SNR 0509-68.7 are confirmed, 
the frequency of ex-companions in known remnants will 
be consistent with the frequency of overluminous SNe 
la compared to SNe la overall. Further work will help 
clarify these events and demonstrate if in fact this is the 
case. 

Additionally, we note that light echoes have identified 
a small handful of galactic remnants as normal SNe la 
- see e..g, Krause et al. (2008). SNR 0509-67.5 is the 
only extragalactic remnant which has a subtype identi¬ 
fication confirmed by light echoes; it has been identified 
as an overluminous supernova remnant (Rest et al. 2008; 
Badenes et al. 2008). The overluminous subtype classifi¬ 
cation of SNR 0509-67.5 presents an interesting challenge 
to the framework presented in this paper. This event ap¬ 
pears to be near the boundary of normal and overlumi¬ 
nous events; of the six best spectral template matches to 
SNR 0509-67.5, three are SN 1991T-like events, and three 
are normal, though the normal events are not as in good 
agreement as the overluminous events. Alternatively, it 
is also possible that double-degenerates may give rise to 
overluminous SNe la alongside single-degenerates. Spec¬ 
tral classification of additional remnants, most crucially 
0505-67.9 and SNR 0509-68.7, if confirmed to harbor ex¬ 
companions, as well as 3C 397 (see below) will yield cru¬ 
cial data confronting theory. 

Lastly, one of the most direct diagnostics of the amount 
of ®®Ni synthesized in a SNe la event is the hard X-ray 
Fe K-a line in young remnants. The hard X-ray portion 
of the spectrum also contains adjacent K-a lines of sta¬ 
ble iron-peak elements including ®®Ni, ®®Mn, and ®^Cr, 
which are produced by electron captures in significant 
abundances only at the high densities of Chandrasekhar- 
mass white dwarfs (Seitenzahl et al. 2013), and are 
a characteristic nucleosynthetic signature of the single¬ 
degenerate channel. Direct detection of the 5.9 keV Mn 
K-shell line produced by the decay of radioactive ®^Fe 
may also be possible in the first several years follow¬ 
ing nearby SN la events(Seitenzahl et al. 2015). After 
the initial submission of this paper, significant levels of 
stable Fe-peak elements in Suzaku X-ray data were an- 
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nounced for SNR 3C 397 (Yamaguchi et al. 2015), yield¬ 
ing the strongest supporting evidence to date of a single¬ 
degenerate SN la. Furthermore, 3C 397 has one of the 
highest measured X-ray luminosities in the Fe K-a line 
of all SNe la remnants, suggestive of an overluminous nu- 
cleosynthetic yield of ^ I.OMq of ^®Ni (Yamaguchi et al. 
2014; Patnaude et al. 2015), though more work is needed 
to compare these recent observations against theoretical 
models and more precisely determine the likely ®®Ni yield 
of this event. A crucial test of the self-consistency of the 
single-degenerate paradigm will be to determine if 3C 397 
harbors an ex-companion. A light echo spectra for 3C 
397 of sufficiently high quality to determine the subtype 
of the event will be a further independent confirmation 
of the luminosity. Theoretical work confronting three- 
dimensional SNe la explosion models, evolved through 
the young SNR remnant phase, with this recent hard 
X-ray spectral data will further shed light on the single¬ 
degenerate explosion mechanism. Additionally, and most 
crucially, further observations and modeling are required 
to detect stable Fe-peak elements in other remnants and 
to quantify their ^®Ni yields, which will provide cru¬ 
cial constraints on the relative contributions of both the 
single-degenerate and double-degenerate channels to the 
overall SNe la rate, and inform optical light curve mod¬ 
els. 

Combining all available observational and theo¬ 
retical evidence, the picture that emerges is that 
Chandrasekhar-mass white dwarfs, long thought to be 
the origin of the standardizable candle, may instead 
be the anomalies. In this view, the stellar progenitors 
of single-degenerate SNe la are primarily the rapidly- 
accreting supersoft X-ray sources, leading to overlumi¬ 
nous SN 1991T-like events, with recurrent novae respon¬ 
sible for the smaller class of SNe la-CSM events. Normal 
SNe la events are very likely to instead originate primar¬ 
ily from a combination of the double-degenerate channel 
and double-detonations. Consequently, the absence of 
evidence for the single-degenerate channel, derived al¬ 
most entirely from the more commonplace normal SNe 
la, is not altogether surprising. All evidence points to 
the conclusion that single-degenerate SNe must be rare 
in comparison to the total SNe la rate, the proverbial 
needle in the haystack. In order to either confirm or 
strongly rule out the single-degenerate channel, one must 
know where in the haystack to look. The class of over- 
luminous SNe la represents the most likely outcome of 
the single-degenerate channel, and future observations 
focused upon overluminous events may help elucidate 
their nature. 
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APPENDIX 

In this appendix, we fill in the physical and mathematical modeling of the birth of a flame bubble during the 
simmering phase of Chandrasekhar-mass white dwarfs. We assume that the flame bubble begins with ignition of a 
spherical bubble with radius Rq. The assumption of sphericity is an excellent one for small flame bubbles less than the 
fire-polishing scale Afp defined in the text, since the action of the flame will rapidly erase any non-spherical distortions. 
In this case, the bubble radius R{t) satisfies the equation : 


R = 


l-bTl 

1 -7l 


Si 


( 1 ) 


where the Atwood number A is defined in the text. Numerical simulations demonstrate that the distribution of offset 
radii occurs well within the central pressure scale height (~ 400 km) of the white dwarf, so the density stratification 
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within the ignition region may be neglected to a good degree of approximation. With the zero time chosen as the 
ignition of the bubble, the bubble radius R is then simply expressed as 

m = Ro+(^\^^s,t ( 2 ) 

where Si is the laminar flame speed. 

The difference between the density of the ash in the bubble and the density of the surrounding fuel will soon cause the 
bubble to buoyantly rise. Additionally, the bubble is subject to the added mass effect, in which the action of displacing 
the surrounding fuel causes upward inertia in the bubble, effectively increasing its buoyancy (Landau & Lifshitz 1959). 
For a spherical bubble, the added mass is precisely 1/2 the mass of the fuel displaced. We can summarize these effects 
with the Morison Equation (Brennen 1995) : 


dt 




iPf 


V 


iPf 


Pa) 9 


(3) 


here t is time, R is the radius of the bubble, pa is the ash density of the bubble, p/ is the fuel density of the bubble, v 
is its rise speed, g is local gravitational acceleration. 

The left-hand side of equation 3 is the time derivative of the inertia. The first term is the inertia of the bubble 
itself, whereas the second is the inertia contributed by the added mass effect, exactly equal to half that of the fuel 
density displaced by the spherical bubble. In contrast, the right hand side gives the buoyancy force. Taking the time 
derivative of inertia yields both an acceleration and a drag term, which arises from the growing mass of the bubble. 
One can easily demonstrate that the net acceleration of the bubble is : 


dv 

dt 


ah - ad 


where the acceleration due to buoyancy, ah, and drag, ad, are : 


ah — 


A 

1 - A 


(4) 

(5) 


ad 


/I -I- A\ Siv 

[i-aJ ~R 


and Ka is a factor taking into account the added-mass effect, and is defined as 


( 6 ) 


Ka = 


1-, 1 1+A 
^ ^ 2 1-A 


(7) 


We note that our equations 1 and 3 are nearly identical to those given by Aspden et al. (2011), but differ in one 
respect. Our model takes the flame bubble growth at the laminar flame speed, which is appropriate to the early- 
growth of the flame bubble. In contrast, Aspden et al. (2011) also include the effect of entrainment, which adds a term 
proportional to the bubble rise speed v to equation 1. Additionally, Aspden et al. (2011) use two free parameters to fit 
the effective buoyancy and drag of their model to numerical simulations, whereas our approach has no free parameters. 
Nonetheless, the two approaches are in good agreement during the initial stages of bubble growth following ignition, 
where the bubble rise speed is less than the laminar flame speed. 

We now derive an equation of motion for the bubble offset radius r(t). Buoyant acceleration is independent of the 
size of the bubble, whereas the drag acceleration depends inversely on bubble radius. Consequently, the motion of 
a small flame bubble is strongly influenced by drag. The ignition points obtained within numerical simulations are 
typically very small (of order km), and are themselves likely limited by numerical resolution. Hence, we are motivated 
to taking the limit of zero initial flame bubble radius, going over to the self-similar limit, with the bubble radius simply 
becoming R{t) = [(1 -|- A)/{I — A)\Sit. 

Writing the gravitational acceleration in the central, nearly-uniform density portion of the WD as g = g^r, the offset 
radius of a self-similar bubble satisfies the following equation : 


'dV 

dt"^ 


-^Ka 


A 

1-A 


gor 


„dr 


( 8 ) 


The equation is singular at t = 0 due to the assumption of zero bubble radius at ignition, but is well-behaved 
provided the bubble begins at rest with respect to the background flow - a reasonable assumption which we adopt. 
The solution to this equation may be found by a power-series expansion in time. Under the assumption of zero initial 
velocity, it is easily demonstrated that the odd terms in this series all identically vanish. To leading order. 



r{t) = ?'o + ^ 


3-A 


(9) 
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This implies an effective acceleration equal to 


A 

3-A 


which is the value we use in the text. 



